Low temperature superfluid stiffness of rf-wave superconductor in a magnetic field 
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The temperature and field dependence of the superfluid density p 3 in the vortex state of a d-wave 
superconductor are calculated using a microscopic model in the quasiclassical approximation. We 



show that at temperatures below T* oc y/H, the linear T dependence of p s crosses over to a T 
dependence differently from the behavior of the effective penetration depth, \~^{T). We point out 
that the expected dependences could be probed by a mutual-inductance technique experiment. 
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PACS numbers: 74.25.Nf, 74.72.-h, 74.25.Jb 

Introduction. While the mechanism of superconduc- 
tivity and unusual nature of the normal state in high- 
temperature superconductors (HTSC) are not yet un- 
derstood, there is a consensus that the zero field super- 
conducting state has a d-wave superconducting energy 
gap, with nodes along the diagonals of the Brillouin zone 
jij. The presence of the nodes results in a density of 
low energy quasiparticle excitations large compared with 
conventional s-wave superconductors even at the tem- 
peratures much smaller than the transition temperature, 
T « T c . Although these excitations are reasonably well 
described by Landau quasiparticles their presence brings 
in a qualitatively new quasiparticle phenomenology not 
encountered in conventional superconductors. Among 
many aspects of this new physics a major role is played 
by these low energy excitations in the mixed (or vortex) 
state [|[ ||, |], D . Since all HTSCs are extreme type-II su- 
perconductors a huge mixed phase extends from the lower 
critical field, if c i ~ 10 — 100 Gauss to the upper critical 
field, H C 2 ~ 100 T. An important property of the vor- 
tex state was pointed out by Volovik |^|, who predicted, 
that in contrast to conventional superconductors, in d- 
wave systems the density of states (DOS) is dominated 
by contribution from excited quasiparticle states rather 
than the bound states associated with vortex cores. It 
was shown that in an applied magnetic field if the ex- 
tended quasiparticles DOS, N{uj = 0, if) ~ \/H rather 
than ~ if as in the conventional case. This result was 
confirmed by specific heat measurements on high qual- 
ity single crystals (?], g|. Subsequently the semiclassical 
treatment of |6j was incorporated into the Green's func- 
tions formalism extended to include the effects of impu- 
rity scattering |)| . Accounting for the impurity scattering 
which violates a simple y/~H dependence has improved the 
agreement between the theory and measurements of the 
electronic specific heat (see Refs. in Q). 

A similar to the Volovik effect, the weak (if < ff c i) 
field response was studied by Yip and Sauls [nj. They 
predicted a direction dependent nonlinear Meissner effect 
also associated with quasiclassical shift of the excitation 
spectrum due to the superflow created by the screening 
currents. Although initially an experimental evidence for 
such an effect was reported, the subsequent experiments 
lfll| did not confirm this effect. 



Another interesting result was obtained by the /iSR 
measurements of temperature and field dependences of 
the effective in-plane penetration depth, A c ff in single 
crystals YBa2Cu30g.95 by Sonier et al. fll^] . At high 
magnetic fields, they observed a flattening of A~g (de- 
fined in these experiments as the width of magnetic field 
distribution) at low temperatures in contrast to the linear 
T behavior expected in a clean d-w&ve superconductor. 
If one assumes is proportional to the superfluid den- 
sity, then such a flattening could, in principle, indicate 
an opening of a secondary gap at the nodes of a rf-wave 
superconducting gap as was already suggested after the 
measurements of thermal conductivity J2|. 

It was argued by Amin et al. |l3| that the simple re- 
lation A^ oc p s is not valid for the penetration depth 
extracted in pSH experiments at finite fields, so that us- 
ing the proper definition p[ of A c ff, which corresponds 
to the /LtSR measured penetration depth, the observed 
behavior of A^ can be explained by a nonlocal London 
model for a d-wave superconductor. The behavior of the 
superfluid stiffness, defined as a second derivative of the 
free energy of the system with respect to a vector poten- 
tial, in this case cannot be addressed by the /iSR. 

Nevertheless the superfluid density in itself could also 
be extracted, for example, from the measurements of the 
low frequency complex sheet conductance as done in two- 
coil mutual-inductance technique measurements on thin 
films Jll], |l6). The superfluid density (stiffness) in this 
case is found from the inductive part a 2 (w) of the con- 
ductivity fTl 



Ps(T) 
m 47re 2 A|(T) 
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where e is the electron charge, m is its mass (or, more 
generally, a mass of the charge carrier) and c is the ve- 
locity of light. It is very important to make a distinction 
between the London penetration depth appearing in 
(|l|) and A e ff deduced from the /iSR measurements |lj, [18| . 

Thus in this work we investigate the influence of 
the magnetic field directly on the superfluid stiffness in 
HTSC using the quasiclassical approximation at ff c i < 
if <C if C 2. Our result is that at low temperatures there is 
also a flattening of p s which is not related to an opening 
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of the gap and is just the result of a Doppler shift of the 
quasiparticle energies in the vortex state. 

Although a quantum mechanical treatment is certainly 
needed for a consistent treatment of the vortex state || , 
the semiclassical approximation always provides a good 
starting point when a new physical property of the sys- 
tem is involved. Furthermore, as argued in [Q for the 
parameter range relevant to the study of HTSC, this ap- 
proximation reproduces the energy spectrum of the near- 
nodal quasiparticles in a vortex state to a high degree of 
accuracy allowing also to include the effect of impurity 
scattering into the analysis. In particular, the thermal 
conductivity was studied using the semiclassical approx- 
imation g m m . 

Superfluid stiffness in the presence of impurities. Su- 
pcrfluid stiffness is given by p2| : 



pT{t,h) 



a/3 
PsO 



(2) 



The normal fluid density, p n in (Q), calculated within 
the "bubble approximation" with dressed fermion prop- 
agators (i.e., with self-energy E due to the scattering on 
impurities included) but neglecting the vertex and Fermi 
liquid corrections, is 



d 2 k 



m J (2tt) 2 
x tr[GA(w,k)G A (w,k)) 



duo tanh 



w v Fa v F p 



2T Airi 
G B (w,k)G fl (w,k)]. 



(3) 



As shown in |22l the vertex corrections can be neglected if 
the impurity scattering potential is isotropic in k-space. 
Likewise the Fermi liquid corrections can be taken into 
account along the lines of Ref.|22 . In (0) Gr^a{uj, k) are 



the retarded and advanced Green's functions 



Gfl,A(w,k) 



(a;± ? r)/ + r 3 g(k)-T 1 A(k) 
{cu±tT) 2 -£ 2 (k) - A 2 (k) 



(4) 



with the dispersion law £(k), the d-wave supercon- 
ducting gap A(k), and the Fermi velocity is = 
d£(k)/<9k| k=kF . The width T = -ImE(w) in (@) is the 
scattering rate due to impurities and other sources, e.g. 
disordered vortex lattice (see below). Following Q we 
assume that around the nodes T is frequency and mo- 
mentum independent, so that the presence of impurities 
is modeled by a widening of <5-likc quasiparticle peaks by 
a Lorentzian 
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with a constant width T (see, however, |2^] where this 
assumption was not confirmed for dopants considered as 
impurity centers). Here E(k) = y / 4 2 (k) + A 2 (k) is the 
quasiparticle dispersion law. 



Employing for the low temperature regime, T <C T c 
the nodal approximation |]22| we obtain that 



PT(T) v f S, 



a l) 
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J, 
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where va = <9A(k)/<9k|k=k F is the gap velocity and 
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with the cutoff energy po (as estimated in |jj po ~ 
1500K). Taking the no impurities limit, r — > one can 
recover from (^|) , (Q) the known linear T dependence |2^] 



p n (T) 2ki2v F 



T, 
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where we set ti = ks = 1. It is essential, however, to take 
into account the presence of impurities which modifies 
the low temperature T dependence of p n from linear to 
— T 2 J2|, Indeed in the limit I » T we get from 
(§), (fthat 



Pn(T) 
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For clean YBCO monocrystals the estimated value 
of the scattering due to impurities r ~ 1 — 2K, but 
in thin films the quadratic temperature dependence of 
p s (T) is observed over a wide temperature range |ll| im- 
plying that the value of To could be bigger, so that for 
our estimates we use Tq = 10K. However, as we will see 
below the obtained dependences on H are in fact even 
stronger for smaller values of Tq . 

Doppler shift. In the semiclassical approach to the 
vortex state the presence of a superflow is accounted 
for by introducing the Doppler shift into the energy 
uj — > uj — e(k, r) ||, [|, where e(k, r) = v s (r) ■ k, and 
v s (r) is the supervelocity field at a position r created by 
all vortices. There are several ways to treat the problem 
depending on the assumptions we made about the vortex 
lattice. 

1. Completely disordered vortex state. Following |2^] 
we assume that the vortex lattice is disordered, which 
is not unreasonable for thin films ^M- Then on the 
length scales large c ompared with the intervortex dis- 
tance a v = <&o/wB (here $o = hc/2e is the flux quan- 
tum and B is the internal field), propagation of quasipar- 
ticles is described by the Green's function averaged over 
the vortex positions {R^}. This averaging can be done 
using the probability density 



P(r?) = (%-v s (r)-k)) {Ri} , 



(10) 
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so that the averaged Green's function 

/>OC 

(G(u,k)) = / d v V{ V )G\Lo- v ,k), (11) 



where G°(w,k) = G(w,k) from Eq. (|) with Y = 0. For 
typical accessible fields 7J c i < i? <C i? C 2 the magnetiza- 
tion due to the vortex lattice is small, and the internal 
field B can be replaced by the applied field, H directed 
along the c-axis. 

Assuming that the vortex positions are random and 
uncorrelated, one would get p7| a Gaussian distribution 
V{ri] = (2nE 2 H )~ 1 / 2 e-" 2 / 2E H with 



W?{«, 8 (^(r)) 



-ft 2 - 2 



(12) 



where we used the expression |v s (r)| = h/2mr (the 
Planck's constant h is restored in these expressions) 
for a supervelocity field created by a single vortex and 
took r = a v . In what follows we will use an estimate 
iS#[K] ~ 30-\/i7frJ obtained for the typical values of 
vf ~ (1-5 — 2.5) x 10 7 cm/s in YBCO. To make an ana- 
lytical calculations possible we follow Ref . J|(| and replace 
the Gaussian probability distribution by a Lorentzian 
V{rj) = ■n- l E H j{if + E%), so that 



(G R (iJ,-k)) = G R (iJ + iE H ,k). 



(13) 



Thus within employed approximation the effect of quasi- 
particle scattering from the supervelocity field v s (r) cre- 
ated by a completely disordered vortex lattice also results 
in a widening of the quasiparticle peaks with a constant 
width Eh- Although this property may not be valid for 
other forms of the distribution functions, we find it useful 
for a physical insight to the problem. 

We note that the used above Lorentzian distribution 
has another property that the average of the product of 
the Green's functions is equal to the product of averaged 
Green's functions, so that the individual Green's func- 
tion can be averaged over the vortex positions. For any 
other type of the distribution function this property is 
not valid and averaging should be done for the product 
of the Green's functions itself (see below). 

Eq. (|l^) may be used to calculate the field induced 
DOS, N(u, H) = 1/(8tt 2 ) / d 2 ktrA(ui, k), so that N(u = 
0,H) oc E h \tl P 1/E 2 h ~ VHlnH /H with H ~ 2500T. 
This is just DOS from |22 with T replaced by Eh- 



Since po ^ Ejj this result agrees with the dependence 
N(u = 0,H) ~ VH obtained by Volovik ||. Further- 
more, as we will see below, the extra In factor is caused 
by the assumption we made about the completely disor- 
dered vortex state. Note also, that the first temperature 
independent term of Eq. (||) may be identified as a DOS 
contribution to p n . 

Assuming also a Matthiessen type rule for the impurity 
and vortex contributions to the lifetime ^7|, we finally 
obtain that the field dependence of the normal fluid den- 
sity is described by Eqs. @, (g) with T = T + E H - 




FIG. 1: The temperature dependences of p 3 (T, H)/p s (0, 0) 
calculated from Eqs. (|), @, (@) and T = F + E H - The 
dotted lines are calculated using the T 2 asymptotic (0) . 



and its asymptotics for T ^> T and T<T are given by 
Eqs. (||) and (|), respectively. These results are presented 
in Fig. [I] (we used for computation some typical values 
of the parameters pgfl , e.g. vf/va — 15), where one can 
clearly see a crossover from T 2 dependence of p s (T) to the 
linear T dependence at T* ~ To + Eh- For a clean sys- 
tem (Fo = 0) and a single vortex averaging this crossover 
of p s {T) was also obtained in p9|] . Thus this behavior of 
p s (T) indeed differs from the effective penetration depth 

Ki ( T ) El which also 

crosses over from the linear T de- 
pendence, but to a T 3 dependence at (for a clean system) 
T* cc y/H. Note also that since the nodal approximation 
was used to arrive at Eqs. (||), (0) one cannot approach 
the region p s (T, H)/p s (0, 0) <C 1, where the dependence 
of the superconducting gap on temperature should be 
taken into account. 

In Fig. H we show the dependence of p s on the applied 
field H . It appears to be not very different from ~ \Jli 
dependence obtained in |30| for the periodic vortex lattice 
at T = 0. The reason for this coincidence becomes more 
clear after we consider a vortex liquid case. 

2. Vortex liquid. To consider a more ordered vortex 
liquid state we adopt a method used in [Q]. In this case 
one introduces a Doppler shift e n (r) = v s (r) • k„ at the 
node k„ (k x = — k 3 and k 2 = — k4 for a d-wave supercon- 
ductor, so that 6i = —63 and e 2 = —£4) to approximate 
the Doppler shift for the entire nodal region. Then if 
we know how to express a physical quantity F in terms 
of the Green function G(o;,k) we can compute its lo- 
cal value F(r) with the local "Doppler shifted" Green's 
function, G(oj, k; e(k, r)) = G(ui — e n (r), k). We then 
approximate the field-dependent measured value F(H) 
by the spatial average which depends on the magnetic 
field H through e„ as § A' 1 J d 2 rF(e 1 {r), e 2 (r)) = 
deide 2 F(ei, e 2 )£(ei, e 2 ), where the first integral is 
taken over the part of a unit cell of the vortex lattice 
(with the area A) in real space where the Doppler shift is 
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tion of a more ordered vortex glass or even lattice state. 

Hence the local quantity J from Eq. (Q) has to be re- 
placed by its local value J(e„) and 



Pn 
m 



vf 

TTVA 



deJ(e)V{e) 



(15) 



One can obtain that the impurity dominated limit (T <C 
|e„| <C To) is described by Eq. (||) with r = To, while in 
the field dominated, T <C Tq <C |e„| regime 
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FIG. 2: The dependence p 3 (T, H)Jp s (0, 0) on the applied field 
calculated from Eqs. (|), @, (@) and T = T + £^ The 
dotted line is obtained using the first term of Eq. (Iq), i.e. 

for p n = vf/(tvva)Eh- 



much smaller than the gap maximum and in the second 
integral we used the distribution function 



£(ei,e 2 ) = A- 1 J d 2 r5(e 1 



v s (r)ki)5(e 2 - v s (r)k 2 



(14) 



where ki and k 2 label two nearest nodes. If £(ei,e 2 ) 
depends on a single variable e 2 + e 2 , then a more sim- 
ple distribution function V(e) — J_ oo dei£(e,ei) may be 
used. Among several candidates for V{e) suggested in 
@ we choose V{e) = l/2E 2 H (e 2 + E 2 H )- 3 ^ 2 which is the 
most convenient for analytical calculations. Using this 
function one can obtain the following field induced DOS: 
N(u = 0,H) oc V# in the strong field, T < |e n | limit 
and N(u! = 0,H) oc H\nHo/H in the impurity domi- 
nated, r S> |e„| regime. This shows that this distri- 
bution function which allows to reproduce the original 
~ VH Volovik's result is indeed better for the descrip- 



Thus irrespectively of the vortex structure we assumed 
the last ~ T 2 term of ([l^) differs only by a factor n/2 
from Eq. (0) taken with T = T Q + Eh ~ Eh for the 
field dominated regime, Eh ^> Tq. It also agrees with 
p9[ , where, however, there is no second term related to 
impurities. The first term of ( |l6|) is already ~ \^H, not 
~ i/ffln-ffo/^ff, exactly as one would expect from the 
DOS calculation and obtains for a single vortex averaging 
p9[ and the vortex lattice as T = |3^] . 

To conclude we have shown using the quasiclassical 
approximation and the simplest form for the impurity 
scattering that the temperature T* ~ Tq of the crossover 
from the linear to the T 2 dependence of the superfluid 
stiffness p s (T, H) increases in the presence of magnetic 
field to the temperature T* ~ To + Eh ■ The dependence 
of p s (T) is different from the effective penetration depth, 
(T) . The experimental measurements of p s (T, H ) are 
necessary to check our predictions against them. 
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